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In contrast with the behavior of the eddies in the open-ocean, the sub-mesoscale 
eddies generated in the constricted Ligurian Basin (NW Mediterranean), are 
unproductive but their combined effect, arranged in a rim-hke fashion, con- 
tributes to the containment of a Productive Ligurian Pool (PLP). Data de- 
rived from MODIS satelhte sensor showed persistent higher chlorophyll con- 
centrations in the centre of the basin, concurrent with high EKE values in 
its surroundings, derived from AVISO altimetry merged products. This sug- 
gested that this 'productive pool' is maintained by the intense (sub)mesoscale 
eddy activity in the rim. Numerical realistic experiments, using a Regional 
Ocean Model System, forced by MERCATOR and by a high-resolution COSMO- 
17 atmospheric model, also showed that most of the sub-mesoscale eddies, 
during 2009 and 2010, are concentrated in the rim surrounding the basin, 
contributing to the formation of a basin-scale cyclonic gyre. We hypothesized 
that the interaction between eddies in the rim might contribute to import 
of nutrients into the pool in two ways: (i) by advection of nutrients from the 
nearby coastal regions into the pool; (ii) by concentrating eddy upwelled nu- 
trients inside the pool; or by a combination thereof. 
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1. Introduction 

The northwest part of the Mediterranean Sea, the Ligurian Sea, is known as an ohg- 

20 otrophic region [e.g. Nezlin et al, 2004; D'Ortenzio and Ribera d'Alcald, 2009], and there- 
fore any supply of nutrients is extremely valuable to sustain the whole trophic chain. 
The biological production of the Ligurian Sea is dominated by a marked seasonal cycle 
[Arnone, 1994] and a significant inter-annual variability [La Violette, 1994; Marty et al, 
2002] . Nevertheless, the system heterogeneity is closely associated with local and regional 

25 hydrodynamical factors, especially those responsible for a high variability of the mixed 
layer depth and for the transport processes [Nezlin et al, 2004; Casella et ai, 2011]. The 
regional hydrodynamics is dominated by a cyclonic circulation system feed by the East 
Corsica Current (ECC) and the West Corsica Current (WCC) which is known as the 
Liguro-Provengal-Catalan (LPC) current (Figure 1). The LPC-current flows along the 

30 costal slopes of Italy, France and Spain, and is affected by instability processes which 
generate (sub)mesoscale eddies, capable of inducing relatively intense shelf-edge flows, 
producing significant dynamical heterogeneity [Millot, 1991]. The presence of mesoscale ed- 
dies in the Western Mediterranean Sea is well documented [Santoleri et ai, 1983; Marullo 
et a/., 1985; Gasparini et ai, 1999; Robinson and Leslie, 2001; Echevin et al, 2003; Casella 

35 et al, 2011]. Nevertheless, the peculiar feature of the Northwest Mediterranean Sea em- 
phasized herein is that, in contrast to many other marginal seas rich in mesoscale eddy 
activity, the nearshore zones are poor in surface chlorophyll [Nezlin et ai, 2004]. The 
LPC-current and the presence of mesoscale eddy activity, is reflected in the high values of 
EKE (Figure 2). In this discussion, we focus on role of a rim of eddies along the Ligurian 
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40 Pronvengal basin, for the maintenance of a cyclonic gyre, and for the containment of the 
Productive Ligurian Pool (PLP, hereafter). 

2. Data and methods 
2.1. Satellite products 

Eddy Kinetic Energy (EKE) was computed from altimetry data. Surface velocities were 

computed from weekly merged products of absolute dynamic topography (ADT), at 1/8° 

resolution on a Mercator projection, distributed by AVISO (www.aviso.oceanobs.com). 

45 ADT maps are obtained by merging measurements from all available altimeter missions 
[Ducet et ai, 2000]. Combining data from different missions significantly improves the 
estimation of mesoscale signals [Pascual et ai, 2006]. We used the processed series which 
considers up to 4 satellites at a given time, thus it has the best possible regional sampling. 
Furthermore, the period of our study (2009-2010) has one of the best altimetric coverage 

50 with a four-satellite constellation (Jason 1, Jason 2, Envisat and Cryostat). ADT is 
obtained adding along track Sea Level Anomaly to the Mean dynamic Topography [Rio 
and Hernandez, 2004]. 

Chlorophyl data was derived from the MODIS-Moderate Resolution Imaging Spectro- 
radiometer flying onboard of the Aqua platform. MODIS data was processed by the 

55 Ocean Biology Processing Group (OBPG) at Goddard Space Flight Center. Merged 
Level-3 chlorophyll products are being created routinely for daily, 8-day, monthly, sea- 
sonal and annual time periods, processed at a 9km spatial resolution. Although seasonal 
and monthly averages were analyzed, only the annual average product is shown, in order 
to emphasize the persistence of the PLP event during 2009 and 2010. 
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2.2. The Ligurian-Provengal Regional Ocean Modeling System 

60 The Regional Ocean Modeling System (ROMS) [Shchepetkin and McWilliams, 
2003, 2005] was configured for the Northwestern part of the Mediterranean Sea. ROMS 
solves the primitive equations based on the Boussinesq approximation. In the Ligurian 
model solutions for 2009 and 2010, the model domain extends from 37.8° N to 44.5° N 
and from 2° E to 16.5° E. The bottom topography is derived from a 30 arc-second resolu- 

65 tion database GEBCO08. In the model version adopted here, there are two open lateral 
boundaries: at 37.8° N and 2° E. We have chosen a horizontal resolution of 1/32°. At 
this resolution, the Rossby radius of deformation (of the order of 5 — 12km in the whole 
Mediterranean and for different seasons, see Grilli and Pinardi [1998] is resolved and 
consequently the model configuration is adequate to simulate mesoscale structures. The 

70 model grid has 35 vertical levels with vertical refinement near the surface, to obtain a 
satisfactory representation of the surface layer and the euphotic zone. At the open lateral 
boundaries, the model is forced with temperature, salinity and velocity fields obtained 
from the MERCATOR product PSY2V3 (www.mercator-ocean.fr). MERCATOR, has an 
horizontal spatial resolution of 1/12°, with daily outputs. 

75 At the sea surface, the regional ocean circulation model was forced with the monthly 
mean climatologies of heat and freshwater fiuxes, derived from the Comprehensive Ocean- 
Atmosphere Data Set, COADS[da Silva A.M. et al, 1994]. For the atmospheric mo- 
mentum, wind-stress was extracted from the Limited Area Model Italy (COSMO-17). 
COSMO-17 [Montani et ai, 2003], is a non-hydrostatic and fully compressible numerical 

80 weather prediction model, which is a regional version of the Lokal Model [Schattler and 
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Doms, 2000], regularly used for operational and research applications. The C0SM0-I7, 
3 hourly solutions, have an horizontal resolution of 1/16°. Validation of the C0SM0-I7 
wind fields has been widely performed [Steppeler et ai, 2003]. As shown in Casella et al. 
[2011], forcing the Ligurian ocean circulation model with high-resolution winds, substan- 
85 tially contributed to the formation of mesoscale and sub-mesoscale eddies, which is an 
important characteristic of the regional dynamics [Gasparini et al, 1999; Robinson and 
Leslie, 2001; Echevin et al, 2003; Nezlin et al, 2004]. 

2.3. ROMS data comparison 

Non-spite the scarce availability of public data, we have compared the two years (2009, 
2010) of ROMS numerical simulations with current meter data, deployed in the Corsica 

90 Channel. Data has been collected by the 'Istituto di Scienze Marine Consiglio Nazionale 
delle Ricerche' (ISMAR-CNR, La Spezia). We considered data from two moorings, located 
at 9.685° E, Longitude and 43.033° N, Latitude, measuring at two different depths (70m, 
125m). From March 2010, the current meters have been replaced by an Acoustic Doppler 
Current Profiler (ADCP), which measured currents from 13m to about 380m, processed 

95 in 20m bins. During the study period, we compared the simulated current in the Corsica 
Channel with current meters and ADCP data. The model equivalent mean current, are in 
good agreement with the measured data, for the same approximate location. The current 
between 70-125m has the same general direction (north), and comparable magnitudes 
(0.4-0.2 ms^^ in winter; 0.1 — 0.2ms^^ during summer), as well as the same seasonal 

100 variability. ROMS also reproduced the general SST and altimetry, monthly and seasonal 
trends, measured with satellite sensors. 
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2.4. Eddy detection algorithm 

In order to detect and track eddies from our numerical solution, the 'Find Okubo- Weiss 

Eddies in Aviso SSH Product tool contained in the MGET tool for ArcGis' [Roberts 

et al, 2010] was adapted to be used with the ROMS output. As defined by Henson and 

105 Thomas [2008], an eddy consists of a region of high vorticity (the core), surrounded by a 
circulation cell of high strain (the ring). Such regions can be detected using the Okubo- 
Weiss parameter [Okuho^ 1970; Weiss, 1991]. The eddy detection algorithm used the 
Okubo- Weiss parameter (Q) computed for ROMS with Romstools [Penven et a/., 2008], 
and filter the candidate eddy cores using various techniques designed to remove small, 

no ephemeral eddies and other unwanted features. These techniques include removing eddies 
that are smaller than a minimum size and tracking eddies through the time-series, and 
removing eddies that are too short in duration. The final analysis classifies eddies as either 
anticyclonic or cyclonic based on the curvature of the sea surface height (anticyclonic is 
concave down, cyclonic is concave up). After a sensitivity analysis for the study region, 

115 a threshold value of Q = ( — 1.5e~^s~^) was chosen, and only eddies with a minimum 
duration of 3 days, and a minimum area of 8 pixels (about 5km radius of the vorticity 
core), were considered. 

3. The Productive Ligurian Pool 

The analysis of yearly mean ocean color data, confirmed the persistent existence of 

productive zone in centre of the Ligurian Sea, during 2009 and 2010. Figure 2 shows the 

120 spatial distribution of the yearly averaged chlorophyll concentration in 2010, (2009 has a 

similar pattern, figure 2-c). The analysis of the concurrent altimetry data for 2010 also 
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suggests that the region surrounding the LPB has the most intense EKE values (Figure 

2-b). The yearly averaged chlorophyll concentration maintains the imprint of the spatial 
distribution during the spring bloom. Figure 2 shows a dynamic region at the rim domi- 
nated by the ECC, the WCC and the LPC-current, with a low chlorophyll concentration, 
whereas the central area (Productive Ligurian Pool), shows a high chlorophyll concentra- 
tion. In fact, the ROMS simulations confirm the presence of the well known cyclonic gyre 
surrounding the LPB, as shown by streams lines for the 2009 and 2010 model solutions 
(Figure 2). The model clearly shows a cyclonic circulation which surrounds the productive 
area, with intense eddy activity in the rim. To note is the inter-annual variability, the 
main patch of chlorophyll is not located in the same location inside the cyclonic gyre, thus 
suggesting the heterogeneity previously discussed in the literature. 

Eddies are located along the path of the ECC, the WCC and the LPC currents, as 
it is shown in figure 2. A highly dynamic rim of eddies can permit the transport of 
nutrients from the eddy-induced vertical mixing and/or advected from a nearby coastal 
regions to be contained in the central, and less turbulent zone (hence the designative name 
of pool). This transport is expected to have a strongest enriching effect when MLD in 
the central zone is shallower and nutrient can remain in the euphotic layer (i.e. during 
spring). Thus it is not surprising that productivity is twice as high during the spring 
bloom, when compared to the yearly averaged values, nevertheless the general pattern is 
maintained. During 2009, 810 eddies were detected in the realistic ROMS solutions of 
the LPB, 316 (40%) were cyclonic and 485 (60%) anticyclonic; whereas during the 2010 
simulation, 935 eddies were counted, 357 (38%) cyclonic and 578 (62%) anticyclonic. The 
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strongest wind-stress peaks were also calculated for the 2010 solution, suggesting that the 

wind plays a role on the generation of (sub)mesoscale features, which in turn seem to 

condition the overall productivity of the region. In fact, the 2010 spring bloom showed 

higher chlorophyll concentration values, compared to 2009. 

Constricted basins and marginal seas in contrast with the open-ocean can not assume 

steady-state conditions, hence the recent developments on strait-marginal seas systems 

theories. Pratt and Spall [2008], used numerical models to test a theory, which explains 

the exchanged occurring between a schematic buoyancy-forced marginal sea and the open 

ocean. The incoming surface layer formed a baroclinic unstable boundary current, that 

circles the marginal sea in a cyclonic manner and feeding heat to the interior, by way of 

mesoscale eddies. As in the Ligurian-Provengal Basin, eddies play a fundamental role in 

the maintenance and in the exchange processes with the marginal sea. Consistent with 

the overall heat and volume balances for the marginal sea, Pratt and Spall [2008] found 

a continuous family of hydraulically controlled states including critical flows, mediated 

by eddies, in a strait between the marginal sea and the open-ocean. Nevertheless, the 

behavior of eddies in a constricted basin and/or marginal sea is not yet fully understood, 

further investigation of the LPB could continue to provide new insights. Kida et al. [2009] 

showed that the eddies accompanying baroclinic instability in the Faroe Bank Channel 

can set up a double-gyre circulation in the upper ocean, an eddy-driven topographic beta 

plume. In regions where baroclinic instability is growing, the momentum flux from the 

overflow into the upper ocean can act as a drag on the overflow causing the overflow to 

descend the slope at a steeper angle than what would arise from bottom friction alone. In 
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contrast, the upper layer of the Mediterranean overflow is likely to be dominated by an 
entrainment-driven topographic beta plume. The difference arises because entrainment 
occurs at a much shallower location for the Mediterranean case and the background po- 
tential vorticity gradient of the upper ocean is expected to be much larger [Kida et al, 
170 2009]. 
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Figure 1. Bathymetry (meters) and general circulation of the Ligurian-Provengal Basin. TC 

- Thyrenian Current; ECC - East Corsica Current; WCC - West Corsica Current; CC - Corsica 
Channel; LPB - Ligurian Provengal Basin; LPCC - Ligurian Provengal Catalan Current; GL - 
Golf of Lion. The circular patterns represent the rim-eddies and the small arrows the nutrient 
entrainment transport process. 
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Figure 2. Top Panel: (a) Yearly averaged chlorophyll {mgm~^), derived from MODIS- 

Aqua; and (b) yearly averaged EKE {m'^s~'^) derived from AVISO altimetry, for 2010. Bottom 
panel: Yearly averaged ROMS streamlines and eddy detected (black spots), representative of 
the averaged circulation for the Ligurian-Provengal Basin for (c) 2009 and for (d) 2010. Mode 
solutions are overlaid onto yearly averaged chlorophyll data 
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